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Abstract 
The   self -consistent  field  equations  are  derived  for  three 
electrons  in  the   state  of  lowest  energy  in  the  field  of  three   fixed 
protons,   in  analogy  with  the  Hartree-F'ock  equations  for  an  atom.     The 
corresponding  problem  for  two  electrons  and  two  protons  is  also   set  up. 
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INTRODUCTORY 

1.  Introduction 

Some  calculations  have  recently  been  made  on  the  recombination 
of  two  hydrogen  atoms  in  the  presence  of  a  third  hydrogen  atom  .     These 
calculations  are  extremely  crude   and  give  values  for  the   three-body  re- 
combination coefficient  of  the  order  of  lO"       cm  /sec,  which  is  perhaps 
1000   times  smaller  than  the  experimental  values.     Even  though  the  values 
obtained  experimentally  are  rather  unreliable,   it  is  hard  to  understsind 
their  magnitude,   as  the  calculated  recombination  coefficients  are  almost 
of  gas-kinetic  order.     But  any  detailed  discussion  is  limited  by  the  fact 
that  the  potential  energy  of  three  hydrogen  atoms  as  a  function  of  their 
separation  is  not  well  known.     The  form  of  the  potential  used  in  the  re- 
cent calculations  was  very  crude,  though  correct  qualitatively;     quantum 

mechanical  calculations  using  atomic   and  molecular  orbitals  have  been 

2 
made    ,  but  they  are  not  in  good  agreement  with  the   scanty  experimental 

data  availaole.     In  fact,  in  calculations  involving  this  potential  it 

has  been  usual  to  employ  an  empirical  expression. 

It  is  of  interest  to  know  this  intersuition  potential  of  three 
hydrogen  atoms  as  precisely   as  possible,  because  it  is  in  effect  basic  to 
all  calculations  of  the  cross-sections  of  processes  involving  the  transfer 
of  translational  and  vibrational  energy   (e.g.,  normal  chemical  reactions); 
it  is  important  also  in  understanding  the  transfer  of  translational  and 
rotational  energy,  although  this  second  aspect  is  not  to  be  developed  here. 

Accordingly,  it  seems  that  it  would  be  worth  while  making  a  self- 
consistent  field  calculation  to  find  the  ground-state  energy  for  three 
electrons  moving  about  three  prx>tons  A,   B,  C,  which  are  taken  to  be  in  a 
straight  line.     This  is  a  major  computational  problem  which  would  require 
a  high-speed  computing  machine . 

It  may  be  useful  to  recall  here  that  a  self -consistent  field  cal- 
culation will  give  the  best  eigenfunctions  and  eigenvalues  that  can  be  ob- 
tained by  any  calculation  using  a  total  wave  function  of  the  form  that  is 
specified:      thus  in  the   case  of  an  atom  with  N  electrons,    the  Hartree  calcu- 
lation gives  the  best  results  for  a  total  wave  function  consisting  of  a  pro- 
duct of  N  one-electron  wave  functions  (i.e.,   without  consideration  of  ex- 
change).    Similarly,   the  Fock  calculation  gives  the  best  results  for  a  total 
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wave  function  consisting  of  N  one-electron  functions  when  the  Fauli  prin- 
ciple is  taken  into  account  (by  using  a  determinantel  total  wave  function)  • 

The  calculation  that  is  contemplated  here  is  an  exact  analogue 
of  the  usual  Fock  calculation  for  atoms:  for  N  electrons  one  has  N  one- 
electron  functions,  determined  by  N  coupled  integro-differential  equations. 
The  exclusion  principle  is  tal-.en  into  account  properly,  and  thus  the  only 
error  (aside  from  spin  and  other  higher  order  effects)  mscy  be  considered 
as  arising  from  the  neglect  of  correlations  of  the  positions  of  different 

electrons,  since  there  are  no  terms  in  r. .  in  the  wave  function,  where  r. , 
'  ij  '       i' 

r.  are  the  coordinates  of  two  electrons.  This  neglect  of  correlation  is 

likely  to  be  fairly  important,  and  the  situation  is  difficult  in  that  the 
"natural"  unit  of  energy  in  calculations  involving  atomic  electrons  is 
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e  /a     »  2.  Ry  »  27.2  ev,  whereas  in  problems  concerning  slow  collisions 

the  much  smaller  mean  gas  kinetic  energy,  kT  -  0 .03  ev  at  300     ,   is  the 
unit  of  importance.     However,  the   calculation  seems  worth  while  in  spite 
of  this,  because  it  is  definitive  insofar  as  it  sets  the  limit  attainable 
by  all  calculations  using  atomic  or  molecular  orbital  functions,   and  would 
be  the  first  cslculation  for  a  simple  moleculrr  system  that  is  strictly  com- 
parable to  the  Fock  calculation  for  atoms. 

In  vievT  of  the  complexity  of  the  problem  for  H^,  the  correspond- 
ing problem  for  H-  has  been  set  up  in  detail  (section  2).     In  the  case  of 
Hp  a  large  number  of  calculations  have  been  made   and  results  in  good  agree- 
ment with  experiment  have  been  obtained.     Neglecting  the   zero-point  energy 
of  the  vibration,  the  binding  energy  is  knoxfl-i  experimentally  to  be  U.73  ev, 
and  this  value  has  been  obtained  by  the  best  variational  calculations   ; 
with  a  neglect  of  correlation  in  the  wave  function  it   seems  that  the  best 
possible  binding  energy  will  be  U.25  -  ii.30  ev.      The  Hp  problem  involves 
all  the    same  procedures  in  physics  and  computation  as  the  problem  for  H^, 
except  that  the  exact  answer  is  known  and  one  has  two  unknown  functions 

In  principle  one  can  get  an  exact  answer  to  the  problem  by  using  a  total 
wave  function  consisting  of  kN  one-electron  wavefunctions   (with  or  with- 
out exchange)   and  choosing  sufficiently  large  integer  k  (in  principle, 
letting  k  tend  to  infinity).     However,  with  kN  one-electron  functions, 
one  would  have  to  solve  kN  coupled  integro-differential  equations,   and 
thus   far  it  has  been  too  complicated  to  consider  k  >  1.     This  point  was 
made  by  Dr.   T.  Y.  Wu. 
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instead  of  three,  and  in  general  the  scale  of  the  work  required  is  ireduced. 
It  would  possibly  be  useful  to  make  the  calculation  for  Hp  for  one  separa- 
tion of  the  protons  merely  to  see  how  rapidly  the  computation  of  the  binding 
energy  converges,  before  starting  on  the  H^  problem;  although  it  would  also 

be  quite  interesting  to  have  the  numerical  value  of  the  binding  energy  of  Hp 
given  by  this  method. 


2.  The 


"2 


Problem 


The  total  ground-state  wave  function,  in  terms  of  one-electron 
wave  functions,  is 

a(l)  .  p(l)  b(l)  .  a(l) 


(2.1)  T  =  K 


a(l)  .  a(l)  b(l)  .  p(l) 
a(2)  .  a(2)  b(2)  •  ^(2) 


a(2)  .  p(2)  b(2)  .  a(2) 


T-jhere  a,  ^  are  one -electron  spin  functions  (corresponding  to  <r  «  +  ^  ti  and 
'2^  respectively),  and  where  a(i),  b(i)  are  one-electron  space  wave  func- 
tions  m  the  coordinate  r.  (i  =  1,2),  In  the  zero-order  approximation  they 
are   (is)  functions  centered  about  the  two  nuclei  A,  B  respectively,   so  that 


a  ( 


0-/4 


)  exp(-Zr. ,):  r. , 
ji  ^         iA  '     lA 


'   1  A' 


In  general. 


a(i)     -  C(i^^)    [exp(-Zr^^)  +  Kl  (r^)exp(-Zr^g)  ]    . 


If  C   =  constant,   then  k,  »  0  corresponds  to  Heitler-London  functions  (atomic 

orbitals),   and  >^  =  1  to  Hund-Mulliken  functions  (molecular  orbitals).     K  is 

a  normalization  constant  detenu.ned  by 

(2.2)  (Y,Y)   '  1; 

we  also  have  the  following  normalization  conditions  on  a(i),  b(i): 


(2.3) 


(a, a) 
1  (b,b) 


1 
1 


which  of  course  mean 

All  wave  functions  used  here  are  real. 


1,  etc. 


-  A  - 

Now,  we  wish  to  set  up  the  variation  principle  for  the  energy, 
with  the  auxiliary  conditions   (2.2),   (2.3)   - 

(2.U)  5    [e-L(Y,T)  -  M(a,a)  -  N(b,b)  ]  =  0   , 

vrtiere 

(2.5)  E  =  (T,HY), 
and 

(2.6)  H  -  X^       H(i)  +  V(l;2);   H(i)   =  -  l^Mrr^  *  —-)',   V(l;2)  =  -^  , 

where  we  have  used  atomic  units,   so  that  e  =  ti  =  m  =  1;     the  unit  of  length 

is  the  Bohr  radius  a     =  1i  /     o  "  0.^29  A,   and  the  unit  of  energy  is 

°  me 

2 
e  /a     =  2Ry  »  27.2  ev.  L,M,N  are  Lagrange  multipliers  which  arise  from  the 

subsidiary  conditions  (2.2),  (2.3).  If  one  substitutes  equations  (2.1), 
(2.2),  (2.3),  (2.5)  and  (2.6)  in  (2.U),  one  obtains  the  basic  wave  equa- 
tions  for  a(i),  b(i)  in  the  following  form 

(2.7)  [H(i)  +  V^^(i)-M  Ja(i)  +    [  (a,b)  |^H(i)-EJ.  +  V^^(i)  +  H^^  ]  b(i)  -  0 

(2.8)  [(a,b)^H(i)-E]+V^^^(i)+  H^^  ]a(i)  +    [H(i)+  V^^(i)  -  nJ   b(i)  -  0 

to  be  solved  with  the  subsidiary  conditions  (2.3).  The  abbreviations  are  the 
conventional  ones: 

"ab  =/H(i)a(i)b(i)  d^f^ 

V^t(i)  ■=/v(i,j)a(j)b(j)d^?^ 

etc. 

Thus,  we  have  to  solve  the  system  (2.l)-(2.8)  for  the  ground  state: 
the  problem  has  cylindrical  symmetry  about  the  axis  of  the  nuclei  AB,  and 
hence  r.  has  just  two  free  variables,  r. ,0.  or  p.,  z.  -  the  azimuth  p  does 
not  enter. 

To  sum  up,  we  wish  to  solve  two  coupled  partial  differential  equa- 
tions, (2.7),  (2.8)  with  two  independent  variables,  p,2  (say)  for  two  un- 
known functions  a,b,  and  two  iinknown  constants  M,N  to  be  determined  from 
the  subsidiary  conditions  (2.3).  The  functions  a,  b  are  finite  at  the  singu- 
larities A,  B  and  everywhere  else;  they  vanish  at  infinity. 
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The  staridard  method  of  solution  is  one  of  iteration:  one  starts 
with  fairly  good  initial  wave  functions  a  ,  b  :  from  these  one  calculates 
V  ,  (i),  H  ,  ,  E°  etc., substitutes  them  in  (2.7),   (2.8),   and  obtains  a  ,  b 


(which  satisfy   (2.3)   as  well).     Then  one  calculates  \r,  (i),   H  ,  ,  E     etc., 

ab  ab 

and  repeats  the  process  until  successive  E-values  agree   sufficiently  well, 

2 
The  unit  of  energy  is  e  /&    ~  2?  ev,   and  one  wishes  to  calculate   the  binding 

energy  D,  where  E  =  -(e  /a  +  D) .  D°,  which  corresponds  to  a°(i)  -  Ce  i^^ 

etc.,  in  the  region  of  l.U  a    is  from  3.7  t)  3.76  ev  (depending  on  the  precise 
value  chosen  for  the   separation  AB).     Probably  the  best  value  of  the  binding 
energy  attainable  by  this  method  will  be  from  U.25  to  U.30  ev, 

THE  H^  PROBLEM 

The   derivation  of  the  wave  equations  for  this  case  is  quite  lengthy, 
and   accordingly  the  full  details  are  relegated  to   sections  7  to  15 j     here 
(in  sections  3  to  6)  we  give  a  brief  discussion  of  the  problem,  referring 
to  the  results  obtained  in  sections  7  to  15. 

3.  Sketch  of  the  Problem 

This  is  the  major  problem,   and  the  one  of  i*eal  physical  interest. 
We  consider  only  the  case  where  the  protons  A,  B,  C  are  in  a  straight  line. 


X 
A 


X 

B 


X 

C 


The  total  wave  function  is 


T  - 


K^T^  •  "2 '2 


+  K^T. 


(3.1) 


a(l)a(l) 
a(2)a(2) 
a(3)a(3) 


b(l)Ml) 
b(2)p(2) 

b(3)p(3) 


c(l)a(l) 
c(2)a(2) 
c(3)a(3) 


a(l)a(l)   b(l)(3(l)   c(l)a(l) 


a(l)p(l) 
a(2)p(2) 
a(3)p(3) 

a(l)a(l) 


b(l)a(l) 
b(2)a(2) 

• 

b(l)a(l) 

■ 

b(3)a(3) 


c(l)a(l) 
c(2)a(2) 
c(3)a(3) 

c(l)p(l) 


thus  finally 

(3.2)  T  "  K  T^       (f\vC  G^)v*'  , 

'^X-l 

where  the  v     are  orthonormal  spin  functions.     F     and  G     are  functions  of 
a(i),  b(i),  c(i)j  in  fact,  G     arises  from  F     by  a  cyclic  permutation  of 
a,b,c.     The  variational  problem  is  set  up  in  exact  analogy  with  the  two- 
electron  case,  and  the  central  problem  is  the  solution  of  the  following 
three  coupled  partial  differential  equations  (again  with  the  two  free 
variables  p,   z  -  we  still  have  axial  symmetry): 

[-6    +  V^^(?)  -  N^  ]  a(?)     +     6^2[A   +V^2^?)  ]b(?)  +  e^3[Z\   +  V^^{^)  ]   c(?)   -  0 

t^^[A  *  \^ir)  ]   a(r)  +   [A  +  W^^{^)  -  }i^]  b(?)  +  6^3 [A   *V^3(?)  ]     c(?)  =  0 

6^3  [a  +  V^3(?)  ]  a(?)  +  623 [A  +  V23(r)  ]  b(?)  +    [a  +  V33(?)  -  N3^  c(?)  -  0 

where  N, ,Np,N^  are  constants  determined  by  the  subsidiary  conditions 

(3.U)  (a,a)  -  (b,b)  -  (c,c)  -  1; 

Za  is  the  two-dijoensional  Laplacian  operator,   and  the  V.  ,(r)  are  made  up  of 
constants  which  cire  calculated  and  potential  functions  which  have   at  most     1/r 
singularities  at  the  positions  A,B,C  of  the  protons. 

We  require  solutions  a,b,c  which  are  finite  everywhere  and  vanish 
at  infinity.     The  method  of  solving  these  equations  is  similar  to  the  Hp  case: 
once  again,   there  is  an  initial  approximation  of  the  form 

a°(i)   -  c'exp  ^-z'r^^  j 

b°(i)  -  c'exp   [-z'r^g  j 

c°(i)  -   

but  here  of  course  no  final  answer  to  the  problem  is  known. 
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U.  The  Variational  Principle 

Equations   (3.1),   (3.2)   are  equivalent  to 

T  -  K  y^       (F^+  K  G^)  Y^ 
X=l 

and  the  total  hamiltonian  is 

(U.l)  H(total)  -     ^     H(i)  +  i    Z     V(i.j)    . 

i=l  i^'j 

The  eigenvalue  we  want  is 

(U.2)  E  -   (T,  H  T  ) 

subject  to  the  condition 

(U.3)  (T,f )  -  1 

and  also  to 

(U.U)  (a,a)  =  (b,b)  =  (c,c)  -  1 

This  is  a  problem  in  the  calculus  of  variations.  We  must  have 
(U.5)  5[e  -  L(T,Y  )  -  N^'(a,a)  -  N^'(b,b)  -  nJ(c,c)J  -  0 

with  respect  to  variations  of  the  constants  K,K;  and  the  functions  a(i),  b(i) 
c(i).  L,N  -N  are  Lagrange  multipliers  determined  by  the  subsidiary  condi- 

tions  (U.3),  (U.U). 

Now,   in  evciluating  (Y,  H  Y )  we  get  only  the  following  terms 

Hq  -   (F^,H(i)F^)    ;   i  -  X                         Jp  -   (F^,H(i)G^)  i  =  X 

Hj^  -  (F^,H(i)F^)    ;  i  /  X                        J^  -  (F^,H(i)G^)  ±  li  X 

V^  -   (F\v(i,j)F^)   i,j   ji  X(i^j)          Up  =   (F\v(i,j)G^)   i,j  ^  X  (i^j) 
Vj^  -   (F^,V(i,j)F'')  i  or  j  -  X(i7«j)  U^^  -  (F,V(i,j)G)     i  or  j   =  X  (i^j)    . 
Detailed  expressions  for  the   functions  H_^ U„  are  given  in  sec- 
tions 9«10  of  the  appendix.     It  is  clear  that  in  terms  of  these  expressions, 
E  .   (Y,  HY)  -  3K^[  (1+  K^TT    )(Hp+  2Hj^+   7^+  2Vj^) 

.2»e(Jp.  2J^.  V^V}' 
where TT  permutes  a,b,c  cyclically,   i.e.,ll    f(a,b,c)   =  f(b,c,a). 
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Carrying  out  the  variations  with  respect  to  K  andvc  gives,  respectively, 
(U.8)  L  -  E 

(U.9)  ^    --    il^^rii^ 

(F,H-E,F) 
The  variations  with  respect  to  a,b,c  involve  l^inctional  differentiations   (of. 
section  7   ).     The  results  of  the  calculation  give  the  wave  equations  which  are 
presented  in  the  next  section. 

$«  The  Wave  Equations 

This  section  is  a  sunuuaiy  of  the  results  of  sections  8-15.     The 
actual  equations  to  be  solved  are 

(5.1)  ^21  i  *  ^22  b  +  Q23  c  -  0 

Si  i  *  *^32  -  *  "^33  -  '  °    ' 
•vrtiere  the  Q .  .  are  given  in  section  1$.     In  atomic  units, 

(5.2)  H  .  H(?^)   -  -  ^       A  i      -  [l/r.^  .  1/r.g  .  1/r^,} 

(5.3)  V=V(r^.)=    l/r,j 
(5.U)                       W  -  H  -  E, 

where 

(U.7)  E  -   (T,  H  T)   -  3K^  1^(1  +K^TT  )(Hp+  2H^+  Vp+  2W^) 

TT  permutes  a,b,£  cyclically,   and  K  is  a  constant  given  by  the   condition  (U.3) 
as  follows 

(5.5)        (¥,T)   =  3K^    [d+K^TT  )(FF)+  2K(FG)j     -1, 

where  the  scalar  scalar  products  (FF)  and  (FG)  are  given  by  equations  (8.3)  and 

(8.5).  The  Q . ,  contain  three  constants  N  ,  N,  ,N  which  are  determined  from  the 
^  ij  a'  b'  c 

subsidiary  conditions  (3.U). 
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6.  Discussion 

6.]  How  general  is  the  wave  function  Y  ? 

This  point  is  discussed  in  section  1  (see  the  footnote  on  p.  2      ). 
T,   as  given  by  the  self-consistent  field  method,  is  the  best  powwible  wave 
function  that  can  be  obtained  as  one  antisjTtimetrized  product  of  one-electron 
wave  functions;     these  functions  are  chosen  so  that   the   space   symmetry  of  the 
problem  as  a  whole  is  taken  into  account  (to   this  end  we  use  Y,    and  Y^  or  the 

factor  >t  ) .     An  exact  result  could  be  obtained  either  by  taking  an  infinite 

set  of  such  Y  functions,   or  by  letting  the  individual  a,b,c  functions  depend 

on  the  coordinates  of  all  the  particles  in  an  arbitrary  way.     The  effect  of 

the  present  approximation  is   seen  in  section  17  for  the  Hp  problem. 

6.2  Comparison  of  the  lengths  of  the  procedures  for  the  H^  and  the 
H-.  problems. 

In  each  problem  H  (of  (2.6)  or  (5.2))  has  to  be  computed  once  only; 
at  each  step  in  the  iteration  a  number  of  functions  and  constants  have  to  be 
computed.     For  H„  these   are: 

(a)  functions:       a,  b  (which  is  essentially  equivalent  to  a  apart  from  a 

space  inversion) 

V         (V,  ,    is  essentially  equivalent) 

!ab 
i .e . ,   3  functions ; 


(b)   constants:       E,      (ab) 

H  ,   (-  H,,  ) 
aa  bb 


ab 
(a^Vb^),    (abVab) 

(a^Vab)     •=   (b^Vab) 


i.e.,  7  constants 


-lo- 


in the  H^  problem  there   are  essential  differences  between  the 

three  functions  a,b,c,  \inless  r     =  r^„   j     in  that  case  a  and  c  are  equi- 
valent apart  from  a  space  rotation. 
Thus  we  have  the  9  functions:     a,  b,  cj     V^^,   V^^,  V^^;  V^^,  V^^,  V^^ 

and  the  6x3+1=  19  constants:     Ej   (ab),   (be),   (ca);  H^g,H^^,H^^; 

H  ,  ,  H,     ,  H     5     (aW),   (b^Vc'^),   (cW)   ',   (a^Vbc),   (b^Vca),   (c^Vab); 

Slu  DC  Co. 

(abVab)    (caVab)    . 

Thus,  for  the  H^  problem  the  calculation  of  every  point  (specified 
by  the  two   distances  of  separation  of  the  protons,  r  _  and  r^„)  will  be  per- 
haps three  times  as  long  as  for  the  corresponding  H^  problem.     To  get  the 
potential  curve  we  require  for  the  H^  problem,  it  would  be  necessary  to  cal- 
culate 10  to   20  points. 

The  possible  points  to  be  calculated  would  range  as  follows: 

r,„  =  1.0  -  2.!?  ;       r^  ■  1.0  -  2,$  independently  (all  in  units 
AB  oL> 

of  a     =  ti  /me    ) ,     also  some  points  should  be  determined  in  the  range 
^AB  "  ■'■•^  -  °*^^   '       ^BC  =  2.5  -  ii  or  7.$. 
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APPENDIX:      DETAILS  FDR  THE  H-   PROBLEM 


7.  Notation 


The  H^  problem  is  exceedingly  complicated  and  since  many  terms 
appear  in  the  calculation  some  simplifications  are  made  in  the  notation, 
as  follows: 

(i)  A  function  having  a  free  variable  r .     is  denoted  by  f  instead 
of  by  f(rj^),  e.g.,   a  for  a(i)   -  air^):     V^^  for  V^^^^^'  -  ^°^  ^^^±>    * 
(ii)  Inner  products  are  written: 
(ab)     ,     H^^   ,   (a^Vbc)  ior ff  ^  r ^^r ^  (ij)a^(i)b( j)c( j) 

(ab)  (abVbc)  iorff^  r^  d^f*^V(ij)a(i)b(i)b(j)c( j)   , 

etc . 
(iii)  We  use  the  following  permutation  operators: 

TT   :     a->  b->  c->  a  ,  etc. 
1  I  ,  :     a->  b  and  b->  a;  c->  c. 

etc. 

(iv)   Ihe  notation  T'  H  denotes:    "replace  each  term  in  the  H-expression 
by  a  corresponding  V-term  as  follows:" 

(jk)H.j->(jK)V.. 


W 


H  -  E   :  W   [(ij)  +  const.(ik)(kj) J   ->  (k^  Vij)  +  const.   (ikVkj)    . 


(v)  Functional  differentiation,   as  used  in  deriving  the  wave  equa- 
tion from  the  variational  expression  (U.5),  is  denoted  as  follows: 

°A  ^(i'^»£)   ^°^   (6/6a)  f(a,b,c)    ,  and  Dg,   D^  for  6/6b   ,   5/6c    . 
D^(a,a)         -  2a 


D^(ab) 

■  k 

°A«ab 

«»  Hb 

Vaa 

=  2Ha 

D^(a^Vbc) 

-  2aV, 
-  be 

D^(abVca)     -  b V^  *  £ V^ 
D^(abVbc)     -  bV^^ 
Etc. 
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8.     The  Wave  Functions  F,  G; 

T  =  K(T^+»c¥2) 

X         n     V  -X     Xuv  .  /   X  jin        -Xu  >s        .      r.        .  . 

Y     =  a'^  a     p     6         }  V.Y  T  )  "6  y       spin  functions 

F^     -     b(l)    [c(2)a(3)  -  c(3)a(2)]   -  a(l)    |b(2)c(3)  -  b(3)c(2)  ] 

(8^a)      F^     =    b(2)    [c(3)a(l)  -  c(l)a(3)  ]  -  a(2)    [b(3)c(l)  -  b(l)c(3)  1 
F^     =     b(3)    [c(l)a(2)   .... 

G^     -     c(l)  [a(2)b(3)  -a(3)b(2)]   -  b(l)    [c(2)a(3)   -  c(3)a(2)  ] 

(8.2b)      G^     «     c(2)  [a(3)b(l)  -   ... 

G^     -     c(3)  [  ... 
i.e., 

g^=TTf^. 

Also  we  have: 

(8.3)  (F,F)  -  (FF)  =  2[l  -  (ca)2  ]  +  2[l  -  (bc)^  ] 

-  U(ab)    [(be),   (ca)   -   (ab)]    . 

(6.U)  (GO)    -  2|l  -  (ab)^]  +  2[l  -  (ca)^]   -  U(bc)    [(ab)(ca)   -  (be)] 

(8.5)  (FG)     -  -2  -2(ab)^  -  2(bc)^  +  2(ab)(bc)(ca)  +  li(ca)^   . 

The  functional  derivatives  needed  later  are: 

(8.6)  D^(FF)  -  Ub    [2(ab)   -  (bc)(ca)  ]  -  Uc    [(ca)  +   (ab)(bc)  J 

(8.7)  Dg(FF)  =  Ua   [2(ab)   -  (bc)(ca)]  -  Uc    [(be)  +   (ab)(ea)  ] 

(8.8)  D^(FF)  -  -  Ua  [(ca)  +  (ab)(bc)  ]  -  Ub   [(be)  +   (ab)(ca)  ] 

(8.9)  D^(FG)  "  b    [.U(ab)  +  2(bc)(ca) J  +  c    [8(ca)  +  2(ab)(bc)  j 

(8.10)  Dg(FG)  =  a   f-U(ab)  +  2(bc)(ea)  ]  +  c   [-U(bc)  +  2(ab)(c^)  ] 

(8.11)  D^(FG)  -  a   g(c^)   +  2(ab)(bc)  ]  +  b    [-U(bc)  +  2(ab)(ca)  ]    . 
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9.  The  Matrix  Elements  H^ V^^: 


(9.1)  Hjj  -   (F^H(i)F^)   6\       (=  F^H(i)Fl)   =    ...    ) 

i^    =  2H^^^    [l  -  (ca)2]  +  2H^^   [l  -  (bc)^  ]   -  UH^^,   [(bc)(ca)  -  (ab)  ] 

K^  =   (F^H(2)F-'-)    =    ... 

(9.2)  =  H^^  +  H^^  +  2H^^    [i  +   (ab)2  J   +  2(ab)H^^-  2H^^    [(ab)(ca)   +   (be)  J 

-  2H^^   [(ab)(bc)  +   (ca)  ] 

Vp  =   (F-'-V(23)F-'-) 
^^•^^  =  2(b^Vc^)  +  2(c^Va^)  -  2(bcVbc)   -  2(caVca)  -  U(ab)    [(bcVca)-(c^Vab)  J 

V„  -  (r'-v(i2)F^) 

a 

=  2(aW)  +  (b^Vc^)  +   (cW)  -  2(ca)(b^Vca)  -  2(bc)(a^Vbc) 

+  2(ab)(c^Vab)  -  2(ca)(abVbc)  -  2(bc)(abVca)  +  2(abVab) 

10  The  Matrix  Elements  J^ U. ,: 

U  N 


J^  -   (FVi)G^)   &)     '   (r4i(i)G^)   -    ... 
(10 


D 
1)  =  -^"bb&  "  ^""^^^l  *  2H^[(bc)(ca)   -  (ab)  ]   +  211^^^    [(ab)(ca)-(bc)  ] 

-  2H       [(ab)(bc)  -  (ca)  3 


J„  =   (F-Sl(2)G-'-)   =    ... 


N 


2„         „      r,    .    /  ^s2-i 


(10.2)  =  -H^Cl  *  (bc)^J  *  (ca)Xb-  »cc&  *  ^^^  >  (^^^"ab'^^^^^b, 

+  H^^[3(ca)  +  2(ab)(bc)  ] 


-  u  - 


Ujj  -  (F-'-V(23)G-^) 


/         X  m     -2(c^Va^)  +  2(caVca)  +  2(ab)    [(bcVca)  -  (c^Vab)  J 

+2(bc)    [(abVca)  -  (a^Vbc)  ]   -2(ca)    []abVbc)  -  (b^Vca)  ] 


(10 .U) 


Ujj  «  (^^¥(12)0^)   -    ... 

=  -(a^Vb^)   -  (b^Vc^)  +  3(ca)(b^Vca)  -  (ab)(c^Vab)  -  (bc)(a^Vbc) 

+   (caVca)  -  (abVab)  -  (bcVbc)  +  2(ca)   (abVbc) 


11.  Derivatives  of  the  Functions  H^ — V.,: 

u        fJ 


(11.1)  D^Hp     -  UaH   [l  -  (bc)^]  +  Ub   ^kab)  -  (bc)(ca)J   +  H^J-Uc[(ca)Hj^^^+(bc)H^  ^ 

DgHp     =Vi^  D^Hjj=  Ua[H  |(ab)-(bc)(ca)j  +  H^J   +  UbH    [l  -  (ca)^  ] 
(11-2)  .  Uc    [(bc)H^^.   (ca)H^  ] 

(11.3)  Dj,Hp    -  -Ua   [(ca)H^^^+  (b°)"ab]   '^^  t^*^^"aa*   ^^''^"ab  ] 
(ll.U)  D^Hjj     -  2aH  +  2b    [(ab)H  +  H^^*  2(ab)H^^-  (bc)H^^-  (°a)\c  ] 

-  2c   [h   |(ca)  +   (ab)(bc)|  +  H^^+   (ab)Hj^^  ] 

(11.5)  Dgli^  .IT   D^H^  -  2a    [jab)H  *  H^*  2(ab)H     -  (ca)H^^-  (bc)H      ]   +  2bH 

ba  -^ 

-  2c    [h  |(bc)  +   (ab)(ca)}  +  H^^  +  (ab)  H^^J 

(11.6)  ''^''n  "  ■'-  ^-   ^^'^^  '   (ai^)(bc)}-H^^  -   (^^)«bc] 

-2b    JHUbc)  +   (ab)(ca)y +H^^+   (ab)H  ^  1   +  UcH    [l  +   (ab)^l 

(11.7)  D^Vp  -  UaV^p  +  Ub   [(ab)V^^  -  (bcVca)  +   (c^Vab)  J  -  Uc   R  ^  +  (ab)V^^l 

(11.8)  DgVjj  "TTd^Vp  -  Ua   [(ab)V^^  -  (bcVca)  +   (c^Vab)  ]+  UbV^^  -Uc    rv^c+(ab)V^l 

(11.9)  D^Vj^  -  -Ua    JV^.   (ab)V^ J   -Ub    jv^.   (ab)V^J.Uc    [v^-\b"2(ab)V^ J 

(UJD)    D^Vj^=     2a    [2Vbb  *  V  '^^^^''bc]"  ^^   ['V   ^^""^^21  '  ^^^^ V^''^^ V^'^^'^l] 

-2c    [(ca)V^^+   ^^'^^^ab"  ^^^^^^^^  *  (abVbc)J 
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D^V.,  =TT      D.V  -  2a    [iv  .+   (ab)V^^  -  (bc)V       -  (ca)V       +   (c^Vab)  1 

B  N         ,        AN       —   \     ab               cc  ca                 be                      J 

+b    |2V       +  V       -  2(ca)V      1   -  2c    r(bc)V       +   (ca)V       +   (a^Vbc)  +  (abVca)    | 

■~     Ia3.              CC                               CHI  "~     I                  3.0                           0.0                                                                I 


ba 
(11.11) 


D(,Vj^  =  -  2a    r(ca)V^^^+   ^^^^^ab  *   (^^^c^)   *   (abVbc)  1 
(11.12)  _2^   R^°)\a  *   ^°^^^ab  *   (^^^c)   +   (abVca)  J   +  2c    p^^+  V^^  +  2(ab)vn 


12,     Derivatives  of  the  Functions  J  —  V'. 


(12.2) 


( 


D^Jp  =  2b    [h  |(bc)(ca)  -  (ab)].  -  H^^  +  (^^^^bc  "  ^^"^^"ca] 

(12.1)  r      r  ^  ,      .        -I 

+2c    Uh  Hab)(bc)   -   (ca)j  +   (ab)H^^  +   (bc)H^+  H^^  +2(ca)H^^  J 

DgJp  »  2a   Fh   ((bc)(ca)  -  (ab)}  -  H^  +   (ca)H^^-  (^°)"ca]   "  ^^^i    [l  -  (ca)^  ] 

+  2c    1^  ^(ab)(ca)   -   (bc)j  -  H^^  +   (ca)H^^j  -  (ab)H^^  j 

12.3)     D^Jj3  =TT     Vd  =  ^^   [-^  f(ab)(bc)-   (ca)}  ^   (bc)H^>   (ab)H^^.  H^^*  2(ca)H^^    1 
ca  '^ 

+  2b   [^  |(ca)(ab)  -  (be)}  -  %^  *   (°a)"ab  '  (^^^"cal 

Vn  =  -  2aH    [l  +   (bc)^]  -  b    [|l(ab)   +  H^+  2(ab)H^^  -  2(bc)H^^] 

(I2.h) 

+  c   f^|'3(ca)   +2(ab)(bc)}  +  3H^g^  +  2(ca)H^^J 

Vn  "  "-   E^^^-  "  "ab"  2(^^^"cc"  2(^^^"cal  +2bH(ca)^ 
-c[(bc)H.H^^^.  2(bc)H^^-2(ab)H^J 

^C^N  '^      Vn  =  -^   i  {^^"^^^   *  2(ab)(bc)]  +  3H^^-   2(ca)H^^] 

(12.6)  ^^  ^  -I  r  2  1 

-b    [^bc)H  +  Hjj^+  2(bc)H^^  -2(ab)H^J  -  2cH    (l+   (ab)     J 
D^Uj^  =  -Ua    [v^^  +(bc)V^^1-2b    [(ab)V^^  +(°a)V  '^^^^^^ca  ^^c^^^^)   '  (bcVca)  J 

+  2c    r(ca)V,  ,    +(ab)V,      +(bc)V  ,    +  2V^^  +   (b^Vca)   -  (abVbc)  1 
—    L         bb  be  ao  ca  _j 

D  n^  =  >2a   r(ab)V       +(ca)V,      -(bc)V       +  (e^Vab)   -  (bcVca)l   +  Ub(ca)V 
,        -.        B  D  —    I  cc  DC  ca.  -J  — 

(12.8)  *-         —  —  —  ^  -, 

-2c     (bc)V^^  +   (ca)V^jj   -  (ab)V^^  ♦(s?Vbc)    -  (abVca)    I 
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(12  9)       °^"^  ^T:a  ""^"^  "  ^-   ["""iS  ^^'^"^  V   (^^^^bc  *   (ca)V^.(b2vca)-(abVbc)  ] 
-  2b    Rbc)?      +   (ca)V  ,    -(ab)V      +   (a^Vbc)  -  (abVca)    | 

-Uc    [v^.   (ab)V^^ 

^A^N  "  -2i   fbb  *   (*^<^>V]   -^      [^V^'^'^^'cc  -2(ca)V^  .   (c^Vab)  ] 

(12.10)  o  -1 

+  c     2V^^  +3(ca)V^^  +2(abVbc)  ♦  3(b'^Vca) 

^B^N  -  -^   P^ab  *   (^^V'^'^^''^  '   (c^Vab)]  .  2b    [3(ca)V^  -V^  -V^~j 
(12.U) 

-  S.  [2Vbc  -^c)V^  -2(ca)V^  .   (a^Vbc)  ] 

Vn  -'n'        Vn  "  -^   r^^ca  *^^°^^\b  ■^2(abVbc)  +3(b^Vca)  1 
ca  L    ^'~  — 

13,     Some  Useful  Functions 

We  shall  need  the  results  of  the  previous  sections  in  the  following 
forms: 

(13.1)  S  =  Hjj  +  2Hj^  +  Vp  +  2V^  -  E(FF) 

(13.2)  T  -   Jj^  +   2Jj^  +  Up  +   2Uj^  -  E(FG) 
and  now  we  shall  need 

(13.3)  D^_^  S,  IT  D^.c  S,      D^.^T   . 

Accordingly,   the   functions  of   (13.3)   are  listed  here.     We  use  the  notation  +1  H 
of  itection  ?,   (iv)  to  denote  the  V-terms  corresponding  to  the  given  H-tenns. 

(13  .U)     D^S  -  Ua  I^H    [2  -  (bc)^]  .  2  [V^.  V^  -   (bc)V^  J  1^ 

*4^|  w[2(ab)«(bc)(ca)3+2Hgjj-t-   2(ab)H^^-(bc)H^^-(ca)H^^     +    "^H    J 
-Uc|w[(ca)  +(ab)(bc)]+H^^+(ca)H^^+(bc)Hgjj+(ab)H^^  +      Th     } 
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(13.5)         DgS 


(13.6)       D^S 


(13.7)TTd^S 


(l3.8)TTDgS 


(13.9)TT  DpS 


(13.10)     D^T 


Ua  jw[2(ab)-(bc)(ca)]+  2H^^+2(ab)H^^-(ca)H^^-(bc)H^^  +    Ve  j 

.  Ub[H    [2-(ca)2].  2[V^*  V^-  (ca)V^]   } 

-  Ucjw   [(bc)  +  (ab)(ca)  ]  +  H^^+   (bc)H^^+(ca)H^+   (ab)H^^  +   T"^    ^ 

-Uajw    [(ca)  +  (ab)(bc)]]   +  H^^+   (ca)H^^^+   (^<^)"ab*   ^^^^^c  *    ^     1 

-  Ub(w    r(bc)+   (ab)(ca)  ]   +  H     +   (bc)H     +   (ca)H  ,  +   (ab)H       *  TH       V 
+  8c^H    [l  *   (ab)2  ]   .V^  .  v^  .  2(ab)V^| 

-Ua  (w   [(ab)   +   (bc)(ca)]  +  H^+   ^^'^^^cc*  ^"^^^^c 

.   (bc)H^^  .       T  H     } 
*Ub[H    [2-  (ca)2]  .  2[V^.  V^  -  (ca)V^]  j 

+kcfi  [2(bc)-  (ca)(ab)  ]  +  2H^^+  ^^"^^^^aa  "  ^"^^^"ab'^^^^a  *  ^H     J 

-i4a  Iw   [(ca)+   (ab)(bc)]   +  H^^t-   (ca)H^^+   (^^)\c*   ^^^^"ab  *  "^^    J 

+Ub  |w    [2(bc)-  (ca)(ab)]  +  2Hj^^  +  2(bc)H^^  -  (ab)H^^  -  (ca)Hgj^  +  Th   } 
.Uc  [h    [2  -  (ab)2]   .  2    [V^H-  V^  .  (ab)V^  J  j 

8a  [h      [l  .   (bc)2  ]  .  V^.  V^  .   2(bc)V^  ^ 
-Ub  /wQ(ab)+   (bc)(ca)  ^+  H^^  +   (ab)H^^+(ca)H^^  +   (^c)H^^  +    ^^H  J 
-Uc  |w   [(ca)  +   (ab)(bc)  3  +  H^a+  (ca)H^^+  (ab)Hjj^+  (bc)H^  +  Th  ^ 

-  Ua  1^   [l  .   (bc)2j   +  V^*  V^*  2(bc)V^^ 

-  2b  |w    [2(ab)   -   (bc)(ca)3  +  2H^^+  2(ab)H^^  -   (ca)H^j^-  (bc)H^^  ♦  Th    J 
+  2c  |w   [U(ca)  +  (ab)(bc)]  +  ^\i^^+  h{ca)ii^^*   (ab)H^^+   (^c)H^  +111     ]■ 
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(13.11)  DgT  -  -2a  |w    [2(ab)   -  (bc)(ca)  J   +  2H^^+  2(ab)H^^-  (ca)H^^ 

-(bc)H^^.        Ve     } 
-Ub  I^H    [l  -  2(ca)2  ]  .  V^*  V^  -  U(ca)V^  ] 
-2c /w    [2(bc)  -  (ca)(ab)  ]   +  2Hj^^+  2(bc)H^^-  (ca)H^j^-   (ab)H^^+ Th    } 

(13.12)  D^T  -  2a    |w   |Ti(ca)  +   (ab)(bc)  j|  +  UH^^+  Wca)H^j^ 

+   (bc)H^+   (ab)H^^  +     T  H      j- 
-2b  |w  (2(bc)   -  (ca)(ab)]   +  2Hj^^+  2(bc)H^^-  (ca)H^^j-  (ab)H^^+   Vr  } 
-Uc  |h    [l  .   (ab)2  ]  *  V^*  V^  .  2(ab)V^'^ 

Ik.   The  Variational  Principle t   Derivation  of  the  Wave  Equations 

From  Section  h  we  get  our  wave  equations  as 
(lii.l)       D^_^  [(Y,Hf)   -  E(?,?)   -  rr(a,a)  -  Nj(b,b)   -  nJ(c,c)]'    -  0. 

Let  us  examine  this  in  more  detail. 

E  =  (Y,HT)  =  K^  ^  (f^vcg^)  y^H(F^+VCG^)  / 

X,ti=l,2,3 

(1U.2)  =  K^  Z  (fV  rtG^)H(F^+ VCG^) 

X=l,2,3 

=  K^  XIx  {(l+><^^n")(FNlF'')   +  2K(F^^G^)| 


because 


X^ 


and 


(a)     Y^'  /  -  5 

(b)    (G^^G^)  -  (TrF^ilTF^)  -TT(FNiF^) 
(fNig^)  -  (gV^)  , 

where  F  ,G     are  real  functions. 

Again,   since 

(U.l)  H  «  Z     H(i)   +  iy  V(ij) 

1=12,3  i7j 
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we  have 

(11;.3) 

Z^(fNiG^)   =3    [V   2J^*  U^  .  2U^  ] 

and  thus   the  wave  equations  (II4..I)  become: 

3K^  D^  ((1  +  )C^T7)  S  +  2)tT]-  -  2N^'  a  =0 

(lii.U)  3K^  Dg  £(l+)C^rr)  S  +  2Kt}  -  2N^'  b  »  0 

3K^  Dg  1*^^*  '^^^^  S  +  2yCT  ]■  -  2nJ  c  -  0    , 
where  we  have  put 

(1U.5) 

T  =   Jp+  2Jj^+  Uj3+  2Uj^  -  E(FG). 
We  simplify  these  by  using  the  following: 
(1U.6)  N^=  NJ/3K^   N^  =  N^^/3K^    ,  N^  =  nJ/3K^  ; 

(IU.7)         (all     f(a,b,c)      ):  D^TT  f(a,b,c)   =  TT  D^f(a,b,c) 


DgTT    f(a,b,c)   =  TT  D^f 
r^TT    f(a,b,c)   =TT  Dgf   ; 


and  then  (lU.U)  become 


(V  ^^TT\)  S  +  2KD^T       -  2N^a    -     0 

(1U.8)  (Dg+IC^TT  V  ^  *  ^^V       -  2N^  ^  "     ° 

(Dp+>C^TTDg)   S  +  2Yt:Dj,T        -  2N^  c  -     0    . 
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The  N  are  Lagrange  multipliers,  detennined  from  the  conditions  (U.U) .  The 
equations  (IU.6)  may  conveniently  be  written  as 

(1U.9)  Q2I  -  *  '^2  -  ■*■  ^3  -   ■  ° 

Si  -  *  *^2  -  *  ^33  -      "  ° 
The  expressions  for  the  Q.  .  are  given  in  section  15. 

15 .  The  Q,  .  Rinctions 


^i 


It  is  convenient  to  simplify  the  notation  again,  as  follows . 
Take  in,n  as  positive  or  negative  numbers  (in  practice  only  ll,+2,+U  enter), 
and  let  i,j,k  run  over  a,b,c.  Then 

P  [^(ij)  +  n(ik)(kj)3 

stands  for: 

(15.1)  W  ^(ij)+  n(ik)(kj)]+  mH^^+  ra(ij)Hj^  +  "^^^^Kj  *   "^^^J^^k 

+  mV_+  m(io)Vj^+  n(ik)V  +  n(kj)V^^+  m(k^Vij)  ♦  n(ikVkj); 

The  last  line  is  equivalent  to  *Ve  • 

Now  the  non=diagonal  Q's  are  very  simple: 

(15.2)  ^2  '   "^1  "  ^^■'■'  ^  ^^   ^^^^^  '  ^^°)(c^)  ] 

-  aK^P[(ab)  +  (bc)(ca)  ] 

(15.3)  Q23  =  Q32  -     "^   [(be)  +  (ab)(ca)  ] 

+U(>C^-  K  )P    [2(bc)   -  (ab)(ca)  ] 

(15.U)         Q^^  "  ^3  "     -^(^  -"^^^P    C(ca)  +  (ab)(bc)  ] 

+U>C  p(U(ca)  +  (ab)(bc)3  • 
The  diagonal  terms  are  the  following: 

(15.5)    0,1  -  u|h  [2-  (bc)^]  .  2  [Vy^.  V^.  (bc)V^]  j 

+  8(vc2-k:)  [h[i  .  (bc)^J  .  Vy^.  V^.  2(bc)V^^ 

-  2N 
a 
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(15.6)        0^2  "  U(l  .k')  [h    [2-(ca)2]  .  2^V^.  V^  -  ^<^^^'^o.]    ] 
-  8K{y&  -  2(ca)']   *  V^*  V^  -  U(ca)V^  j- 


-2\ 
(15.7)       Q33  =  8(l->C)|ji   [1+  (ab)^]  +  V^^  +  V^^+  2(ab)v^^  \ 

*  U  >c'   I^H   [2.   (ab)2]  .   2    [V^.  V^  -   (ab)V^  Jl 

-  2N      . 
c 

16.      the  H„  Problem  in  the  Present  Simplified  Notation 

The  total  wave  function  is 
(2.1)  Y  «  K[a(l)b(2)  +  a(2)b(l)  J 

and  the  variation  principle  is 

(2^U')  D^^g  l(T,HY)  -  E(7,T)   -  Nj(a,a)   -  Nj(b,b)\     -  0 

where 

(2.6)  H  -  ^    H(i)  +   V(l,2) 

i=l,2 

and  the  N '    are  Lagrange  multipliers  determined  by  the   subsidiary  conditions 
(2.3)  (a,a)  =  (b,b)  =  1. 

We  have  the  following  results: 

(16.1)  (Y,T)   -  2K^    [1  +   (ab)^] 

(16.2)  Z     (Y^HY  )   -  2k2    [h^^.  H^^.  2(ab)Hj 
i=lj2  -* 

(16.3)  (Y,VY)     =     2K2[(a^Vb^)  +   (abVab)  ] 
and  thus,   if  we  write   (2,U')   as 

(16.U)  %  f:    *     %2  -    ■     ° 


Then 
(16.5) 
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Qll  •  «  ^  \b  -  "a 


Q22  -  H  +  V_.  -  W^ 


aa 


where  W  »  H  -  E,   and  E  is  given  by  equation  (2.5)  as  before. 


17.  Data  Available  on  the  Problem 

2,3,U 
This  is  a  summary  of  the   "classical"  resvilts  .     They  are 

given  in  the  original  units,  and  the  numerical  values  should  be  checked  back 
against  the  original  articles.     The  following  notation  is  used: 
H-L  for  Heitler-London: 

3/2 
(17.1)         a  »  -i-  (i.)         exp[-   |?  -  ?^|   /a\  ;  b~exp(-   |?  -  ?g|   /a  \ 


Wang: 
(17.2) 


(17.3) 


.,    3/2 


—  (- — )        exp  4-  Z'  |r  -  r    I   /a   y  ;     b,  c,   similar  (rel.  to  B,C) 
T"    %  L  A        oj         -    - 


^/tT 


Z'   is  a  variational  parameter.     We  have  the  following  units. 
0.529  A. 


a  =  n  /me 


e  /a  -  27.2  ev  =  627  Kcal  /mole. 


17 J.  Hp  problem. 


^AB^^^ 

D  (ev)  -  the  bind 

H-L 

0.80 

3.2 

Wang 

0.76 

3.76  (Z-  -  1.166) 

probable  optimum  of  present  method 

U.27-U.30 

best  calculation 

0.7U 

U.722  1  0.013 

empirical  value 

0.7395 

U.73  +  O.OU 
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17.3  H^ 

Problem  . 

Nuclear  ! 

separations 

2'W% 

Z 

'^bg/% 

2.0 

2.5 

1.5 

2.5 

2.5 

2.5 

2.0 

2.0 

1.5 

1.5 

Binding,  energy  in  Kcal  /mole 
H-L     (Z'-l) 


53.1*1; 
5U.70 
hh.06 
53.11 
IU.99 


Want 

57.23(Z'-1.066) 
57. 56(Z' =1.0565) 
U7.3U(Z' =1.061) 
56.16(Z'»1.059) 
lU.99(Z'-l.00l) 


Hj  ,  for  comparison 


'^AB  "  ^-^5  ^o 


72.18 


^AB-^-^8  % 


86.9U(Z'-1.666) 


In  all  cases  the  binding  energy  is  zero     for  3  infinitely  separated  H-atoms. 
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